Collective Shape-Phases of Interacting Fermion Systems 



Yu-xin Liu,^'^ Zhan-feng Hou,^ Yu Zhang, -'^ and Haiqing Wei'^ 

^Department of Physics and the CMOE Key Laboratory of Heavy Ion Physics, Peking University, Beijing 100871, China 
^Center of Theoretical Nuclear Physics, National Laboratory of Heavy Ion Accelerator, Lanzhou 730000, China 
'^School of Information Science and Engineering, Lanzhou University, Lanzhou 730000, China 

(Dated: February 6, 2008) 

A microscopic theory is presented for identifying shape-phase structures and transitions in inter- 
acting fermion systems. The method provides a microscopic description for collective shape-phases, 
and reveals detailed dependence of such shape-phases on microscopic interaction strengths. The 
theory is generally applicable to fermion systems such as nuclei, quarks, and in particular trapped 
cold atoms, where shape-phases may be observed and investigated in a controlled manner. 
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Collective motion is common in fermion systems 
such as nuclei and molecules. Nuclei have been found 
for a long time to possess interesting geometric shapes, 
such as a vibratingspheroid, an ellipsoid, and exotically 
deformed shapes [H, Pairing induced deformations 
(PIDs) have recently been observed in trapped fermion 
systems of cold atoms [1, Mean- field theories have 
been used to describe the properties of such fermion 
systems as many-electron atoms, nuclei, quarks, and 
trapped atomic fermions @, @, 0] ■ A long-standing ques- 
tion is to identify the shape phase structures and tran- 
sitions in the general framework of mean-field theory in 
fermion space. Only recently there have been studies on 
nuclear shape^hase transitions [1, Q and critical point 
symmetries [lO| in the framework of shell model hav- 
ing either monopole-pairing and quadrupole-quadrupole 
(QQ) interactions or monopole-pairing and quadrupole- 
pairing and on the shape coexistence and evolution in 
superheavy nuclei in the energy density functional the- 
ory involving monopole-pairing term [2] . But the overall 
structure of nuclear shape-phases is still unclear, espe- 
cially in a more general mean-field model that incorpo- 
rates not only monopole-pairing and QQ interactions but 
also one-body terms and quadrupole-pairing. In similar 
systems of trapped atoms, there is also a lack of theo- 
retical understanding of PIDs observed in recent experi- 
ments 

Starting from a microscopic theory of interacting 
fermion systems, this paper presents a systematic method 
for identifying shape-phase structures and transitions in 
such systems, so to provide a unified microscopic de- 
scription for collective shape-phases of both nuclei and 
trapped cold atoms, as well as other similar systems. 
From a microscopic viewpoint, our method recovers the 
rich structures and transitions of the shape-phases (dy- 
namical symmetries) that have been experimentally ob- 
served and accounted for by the interacting boson model 
(IBM) [ul, [l2|. More importantly, our method reveals 
detailed dependence of such shape-phases on the inter- 
action strengths in a mean-field Hamiltonian. The same 
theory may be applied to other systems such as trapped 



cold atoms with tunable interaction parameters, there- 
fore points to the possibility of observing and controlling 
shape-phases and transitions therein. 

As an approximation of the shell model, a successful 
mean-field theory of nuclear structure, the IBM has been 
proven highly successful in describing low-lying collective 
features of medium and heavy mass nuclei [ij, In 
particular, it has well accounted for nuclear shape-phase 



structures and transitions [Tl|, [H, Q [H, [H, [IE 0, [3, 

|20|). The success can be attributed to the incorporation 
of collective degrees of freedom and the simplification 
of calculations by a mapping from fermions to bosons. 
Three dynamical symmetries, U(5), SU(3) and 0(6), are 
naturally incorporated into the IBM, which correspond 
respectively to shape-phases of a spheroid, axially prolate 
rotor and 7-soft rotors U, 13]. There is also SU(3) sym- 
metry corresponding to an axially oblate rotor phase[12|. 
A widely used mean-field Hamiltonian in fermion space 
can be written as 1211 
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where Sj is the single-particle energy, : : denotes the nor- 
mal product of fermion operators, go, 32, k are strengths 
of monopole-pairing, quadrupole-pairing, and QQ inter- 
action, respectively, djm = (~l)"'~™ai-m and = 
r^Y2^. Besides the energy spectrum, a particularly inter- 
ested property is the electric quadrupole transition rate, 
B{E2,U - Lf) - ^^{Lf\\f{E2)\\L,)\ with f{E2) 

being Q2 multiplied by an effective charge. 
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We take a nuclear system as a numerical example, 
although the same theory and numerical procedure can 
be applied to other fermion systems equally well. To sim- 
plify the calculations, the method of Dyson- type boson 
mapping (2^.[23j is employed with the SD pair truncation 
keeping only pairs with angular momentum J = and 
J — 2 Then the Hamiltonian takes an IBM form 

Hb{S'' ,S, D\D), with the S, D pairs being regarded as 
s, d bosons, respectively. Protons and neutrons are not 
differentiated, so the model is actually corresponding to 
the IBM-1. The single-particle wave functions are chosen 
to be harmonic oscillator's with oscillation constant &^ = 
1.0^1/3/^2, A = 130. For the 50-82 shells, the active 
single-particle orbits are 2^5/2,137/2,351/2,1/111/2,2^3/2 
with energies 0,0.8,1.3,2^2.8 MeV, respectively, simi- 
lar to those used in Ref. [25| . The final truncated Hamil- 
tonian is diagonalized in the U(5)-symmetric bases with 
total number of nucleon pairs being set to N — 10. 

To identify the shape-phases, it is helpful to examine 
the correspondence between the interaction strengthes 
in the microscopic model and the dynamical symmetries 
in the IBM. Quantities of interest are the energy ratios 



i?42 = "l^, Rq2 = -gf^, ^02 = "1^, and R 
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0, and ratios of the electric quadrupole 
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transition rates S42 = B\ma\Zl\y -^64 = b(E2;2i^g^) ^ 

o B(E2-fi-2^2i) , o B{E2;22~*2i) ■ , 

^02 = B(i-;2;2i^0i) ' and B22 - B(i;2;2,-.0i) ' ^^ich are 
known to be able to characterize the low-lying energy 
spectrum well. Table H] lists values of these quantities in 
the dynamical symmetries of IBM-1 (with total boson 
number N — 10), which are to be compared with those 
obtained in the microscopic model. 



TABLE I: Values of interested quantities of a 10-boson system 
in the IBM. Those marked with a star depend on additional 
parameters in the Hamiltonian (as given in Refs. jiS, 19]). 
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We first look at the effect of monopole-pairing with 
parameters 50 G [0,0.50] and g2 = k = 0. Figs. [1] and [2] 
show the dependence of the energy levels and the B(E2) 
ratios on go. The degenerate U(5) levels, such as the 
22 — 4i doublet and the O3 — 3i — 42 — 61 quartet, are well 
reproduced at go > 0.12, and the energy and B{E2) ra- 
tios are U(5)-symmetrically valued. So a large go favors 
a spherical phase. The critical value of 50 ~ 0-12 agrees 
well with the enipirical value go ~ 20 /A with A = 130 
assigned in Ref. [26|. In addition, as go is very small, the 
calculation gives approximately the feature of the SU(3) 
or SU(3) symmetry. It indicates that a fermion system 
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FIG. 1: The dependence of low- lying levels on go when (72 = 
A; = 0. 
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FIG. 2: The dependence of -B42, -B64, B02, and ,622 on go 
when g2 — k — 0. 



with very weak monopole-pairing may appear in an axi 
ally deformed shape. 

Next we examine the effects of the QQ interaction, 
which is known to coincide with Elliott's model 
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where the Q2 is a quadrupole tensor in the SU(3) symme- 
try. However, a system may not be limited to the SU(3) 
phase in a more realistic model with one-body terms. 
So we did calculations with parameters k G [0, 0.50] and 
9o = 92 = 0. The energy levels and the B{E2) ratios 
as functions of k are shown in Figs. [3] and [4] respec- 
tively, where k is limited to below 0.20 in part (a) of 
Fig. [3] for a better view of the energy levels. There is 
clearly a special point k « 0.10, at which all the quan- 
tities defined in Table [J are approximately symmetric. 
More specifically, this point corresponds to an 0(6) phase 
with i?42 ~ R22 ~ 2.5, R02 ~ R02 « 4.5, B42 « 1.38, 
Be4 fa 1.52, B02 ~ 0, and ^42 ~ 1.38, in excellent agree- 
ment with the corresponding values in Table HI In ad- 
dition, the degeneracy of the 22 — 4i doublet and the 
O2 — 3i— 42 — 61 quartet is clearly seen in part (a) of Fig. [3] 
On the other hand, at the two extremes of fc < 0.005 and 
k > 0.2, the axially symmetric rotation phase prevails 
and the calculated energies and B{E2) ratios well repro- 
duce those of the SU(3) symmetry in Table IT] 

It is known that an axially symmetric rotor could 
have either a prolate or an oblate shape, corresponding 
to the SU(3) or SU(3) phase in IBM-1. Since the pa- 
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FIG. 3: The calculated energy levels as functions of k when FIG. 5: The dependence of low- lying levels on 32 when go 
go ~ 92 — 0. 0.15 and fc = 0. 
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FIG. 4: The calculated B{E2) ratios £42, Be,4, B02, and B22 FIG. 6: The dependence of B42, 564, B02, and B22 on (72 
as functions of k when go = 52 = 0. when go = 0.15 and fc = 0. 



rameters in Tabic U are the same for both phases, we 
have to implement other quantities to distinguish them. 
One parameter serving this purpose is the quadrupole 
moment (9(2^) which is negative-, zero-, and positive- 
valued in the SU(3), 0(6), and SU(3) phases respectively 
[lit . As the QQ interaction strength k increases from 
below to above 0.1, the calculated Q{2^) (with the ef- 
fective charge normalized to 1) decreases from a positive 
value to zero and then to negative, which indicates an 
evolution of symmetry SU(3) — > 0(6) — * SU(3), namely, 
transitions from oblate to prolate shape with a critical 
point of 7-soft rotation in between. This agrees excel- 
lently with the variation behavior against the parameter 
X in Ref. [l^. In addition, reexamining the effect of the 
go, we know that, as go is very small, the Q{2^) is pos- 
itive. It manifests that the fermion system with weak 
monopole-pairing takes an oblate shape, i.e., it appears 
as a plateau. Such a result reproduces the recentlyob- 
served density distribution of trapped cold atomsp, 0] 
qualitatively. 

Finally, we study how quadrupole-pairing affects the 
shape-phase structure, using parameters 32 G [0,0.15], 
go — 0.15, k — 0. Note that 172 < go holds true for 
most realistic depictions of nuclei [l3|. The calculated 
dependence of low-lying levels and B{E2) ratios on the 
quadrupole-pairing strength is shown in Figs. [H] and [H] 
respectively. For g2 £ [0,0.04], the system is evidently 



in a vibration phase. For 52 £ [0.07,0.10], it is in an 
axially symmetric rotation phase, with the correspond- 
ing parameters in Table [J well reproduced. Further- 
more, (5(2]^) is negative in this region. It shows that 
the (?2 G [0.07, 0.10] generates an axially prolate deforma- 
tion (in SU(3) symmetry). Besides, the 52 e [0.04,0.07] 
induces a U(5)-SU(3) transition, and the 50 > 0.10 sup- 
presses collective motion. 

The above analyses may be summarized in Table Ull 
for the correspondence between regions of interaction 
strengths in our microscopic model and the dynamical 
symmetries {i.e., shape-phases) in the IBM. The Table 
indicates that the occurrence of an axially oblate rotor 
should be rare in nuclei, as it is inaccessible to most of 
the interactions. But it may be common in trapped 
cold atomic systems since the monopole-pairing there 
is very week. It may be noted that, calculations using 
different configuration parameters may result in differ- 
ent regions and critical values of interaction strengths 
corresponding to various shape-phases, but the general 
correspondence between regions of interaction strengths 
and shape-phases and the overall shape-phase structure 
should remain for the characteristic features of an inter- 
acting fermion system. 

In conclusion, we have presented a microscopic 
method of interacting fermion systems and studied the 
dependence of their shape-phases on the strengthes of 
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TABLE II: The correspondence between microscopic param- 
eter settings and shape-phases (dynamical symmetries). 



U(5) SU(3) 0(6) SU{3) 

go (with g2 = k — 0) > 0.12 — — very small 

(with 30 = 0.15, fc = 0) [0,0.04] [0.07,0.10] — — 

k (with ffo = ff2 = 0) — > 0.20 0.10 ^ 0.0 



basic interactions. Our study yields a complete shape- 
phase structure of such systems by providing an extended 
Casten triangle exhibiting modes of collective vibrations, 
axially prolate, oblate, and 7-soft rotations, correspond- 
ing to different regions of monopole-pairing, quadrupole- 
pairing, and QQ interaction strengthes. Specifically, 
strong or weak monopole-pairing leads to a vibration 
or an axially oblate deformation, respectively. While 
quadrupole-pairing and QQ interactions also induce rota- 
tions (i.e., deformations). In detail, a continuous increase 
of the QQ strength sees oblate, 7— soft, and prolate ro- 
tors subsequently. With a fixed strength of vibration- 
corresponding monopole-pairing, a suitable quadrupole- 
pairing strength may induce a transition from spheri- 
cal to axially prolate elliptical phases. But an exceed- 
ingly large quadrupole-pairing strength suppresses col- 
lective motions. Other than nuclei, the theory should be 
applicable to general fermion systems such as trapped 
atomic fermions, where the interaction strengths may be 
conveniently tunable. In particular, our calculations in- 
dicate that an axially deformed shape may be induced 
in trapped atomic systems by either a weak monopole- 
pairing or a suitable quadrupole-pairing interaction, then 
the recently observed PIDs and nonuniform density dis- 
tribution in trapped atomic systems [3, [Ij could be un- 
derstood as the effect of the deformation. Besides, the 
quadrupole-pairing could also be a trigger to a normal to 
LOFF phase transition in fermion systems because the 
-D-pair is just the polarized (i.e., inhomogeneous) pair. 
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